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Abstract 

The QCD sum rule method is apphed to derive a formula for the p, u, meson-nucleon 
spin-isospin-averaged scattering lengths ap^^^^j,. We found that the crucial matrix elements 
are {qj^D^q)]\r (g = u,d) (twist-2 nucleon matrix element) for a^^^ and ms{ss)N- for a^, and 
obtained = 0.14 ± 0.07 fm, = 0.11 ± 0.06 fm and = 0.035 ± 0.020 fm. These 
small numbers originate from a common factor l/{mN + mp^i_j^^). Our result suggests a slight 
increase (< 60 MeV for p, u, and < 15 MeV for 0) of the effective mass of these vector 
mesons in the nuclear matter (in the dilute nucleon gas approximation). The origin of the 
discrepancy with the previous study was clarified. 
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The operator product expansion (OPE) provides us with a convenient tool to decompose a 
variety of correlation functions into the perturbatively calculable c-number coefficients and 
the nonperturbative matrix elements. In its application to the QCD sum rules (QSR) [0, 
H, the OPE supplies an expression for the resonance parameters in terms of the vacuum 
condensates representing the nonperturbative dynamics in the correlators. In the application 
to the deep inelastic scattering (DIS) the OPE isolates the quark-gluon distribution 
functions of the target from the short distance cross sections. In this paper, we investigate 
the vector meson (p, uo, 0)-nucleon scattering lengths utilizing these two aspects of the OPE. 
These scattering lengths can be measured through the photo-production of these vector 
mesons. Furthermore, they determine the mass shift of the vector mesons in the dilute 
nuclear medium. This will be discussed in the last part of this paper. A similar idea was 
recently presented for the nucleon-nucleon scattering length in ref. 0. 

We start our discussion with the forward scattering amplitude of the vector current 
{y = p, uj, (f)) off the nucleon target with the four momentum p = {p^,p) and the polarization 



where q = {ui^q) is the four-momentum carried by and the nucleon state is normalized 
covariantly as {p\p') = (27r)^2p°(5(p — p'). We set p = {mN,0) throughout this work and 
suppress the explicit dependence on p and s. The vector current is defined as J^''^{x) = 
{l/2){u'y^u{x) ^ d'y^d{x)), J^{x) = S7^s(x). Near the pole position of the vector meson, T^i, 
can be associated with the T-matrix for the forward V — N helicity amplitude, %iH,h'H'{i^, q) 
{h{h') and H{H') are the helicities of the initial (final) vector meson and the initial (final) 
nucleon, respectively, and they take the values of h, h' = ±1, and H, H' = ±1/2) as 

ef\q)T,Au^, q)el''>iq) ^ " (^2 Wh^c, q), (2) 

where we introduced the coupling /y and the mass my of the vector meson V by the 
relation (0| V"*-'*-'(g)) = fv^v^^lf^{(l) with the polarization vector t^lf^{q) normalized as 
Spoi. ^f^\Q.)^'^u'' {q) = ~9fiv + 1^i(lu/(f- As is well known in DIS, T^^ can be decomposed into 
the four scalar components respecting the current conservation and the invariance under 
parity and time-reversal. (Two of them correspond to the spin-averaged structure functions 
Wi and W2, and the other two to the spin-dependent ones Gi and G2-) Correspondingly, 
there are 4 independent helicity amplitudes for the vector current-nucleon scattering; T^i ^1, 
T^^^i:zi, '^i.oij — ,oi' ^^^^ being obtained by time-reversal and parity from these 

four. Since the information on Gi and G2 is still lacking, we shall focus on the combination 
T = Ti + {1 — {p ■ qY/m%q^)T2 (ImTj ~ Wi, i = 1,2), which projects the V — N spin- 
averaged T-matrix, T{uj, q). In the low energy limit (g — > 0), T is reduced to the V — N 
spin-averaged scattering length ay = (l/3)(ai/2 + 203/2) (ai/2 and 03/2 are the scattering 
lengths in the spin-1/2 and 3/2 channels, respectively) as T(my, 0) = 247r(m7v + my)ay p. 
A useful quantity for the dispersion analysis is the retarded correlation function defined as 



T^Au,q)=i[ dSe''^-e{x''){ps\ [j;(x), jj(0)] \ps) = - T dn^^^^^S^, (3) 
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which is analytic in the upper half cj-plane with a fixed q. Noting the crossing symmetry, 
the V — N scattering contribution to the spin-averaged spectral function at g = can be 
written as 



-Im 



TT 



e{uj)-\m [T{uj, 0)] - e{-uj)-\m [T{uj, 0)] 



TT 



TT 



2ATTf^m\,{mN + mv)av e{uo)5'{uj^ - ml) ~ e{-uj)5'{uj'^ - ml) +(S.P.), (4) 



where 6'(x 



is the first derivative of the 5-function (double pole term) and (S. P.) denotes the 
simple pole term representing the off-shell energy dependence of the T-matrix. Equation (^) 
can also be derived starting from the spectral representation. Using this form of the spectral 
function in eq. (|]), and noting that the retarded correlation function becomes identical 



to the causal correlation function T^^ in the deep Euclidean region u 
gets 



oo, one 



-2ATcfym'l{mN + mv)av- 



imt 



+ Q 



2^2 



+ 7^os(g') + 7^,(Q' 



(5) 



where we have used the fact that T becomes a function of u"^ in this limit. In eq. @, we 
assumed the spectral function is saturated by the V — N scattering (with its off-shell effect) 
and the "continuum" contribution]^: TZosiQ"^) denotes the simple pole term corresponding 
to the off-shell part of the V — N T-matrix 



1 1 {ml + Q^)) and 7^c((5^) stands for the 



"continuum" contribution with its threshold S'q (~ l/(5'o + Q^)) 0. The sum of the residues 
of TZosiQ"^) and TZdQ^) is constrained by the l/Q^-term in the OPE side of the correlator. 
(See eq. (|D below.) 

We now proceed to the OPE side of T(uj'^ = — Q^) (l.h.s. of eq. (^)). Unlike in DIS, our 
OPE is the short distance expansion and hence all the operators with the same dimension 
contribute in the same order with respect to (= — at g = 0). The complete OPE 
expression for T(Q^) including the operators up to dimension=6 is given in ref. in the 
context of the finite temperature QSR. For the p and u mesons, it reads from eq. (2.13) of 
PI as (— for p and + for u) 



2m, 



2TTa^ 



m 



2Q 



B 



IT + + 



-B- 



(6) 



where denotes the spin- averaged nucleon matrix element, and and are the 
scalar four-quark operators familiar in the QSR for the p and u mesons; = (m7^75A"m =f 
J7^75A'^c?)2, Q+ ^ {u-i^\''u + d^^X''d)Y.Y''qit.>^''q- In eq. (|), Al^" = Al + {n = 
2, 4) are related to the twist-2 operators and are given as the n-th moment of the parton 
distribution function (g = s); {ST {q'j^^D^^ ■ ■ ■ D ^^q{p))) n = (-i)""M^(/i)(p^^ ■ ■ • 
p^^ — traces), A'^{p) = 2 Jq dxx^'^ {q{x, p) + (— l)"g(x, /i)) with the renormalization scale 
p. Bi {i = 1=1=, 2, 3) are associated with the twist-4 matrix elements as {0^^i^{p))n = 



{p^p, - m%g^,/A)Bi{p) with O^J = {gyA)ST [{u^^-f^X'^u T d^t.l5)^''d){p 



Q2 
^ fjLV 
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where the color matrix A° is normahzed as Tr(A'^A'') = 25°''^ and the symbol ST makes the 
operators symmetric and traceless with respect to the Lorentz indices. 

To get an expression for the V — N scattering length, we first make a Borel transform 
of eqs. (^) and (|^) with respect to Q^, and then eliminate the unknown parameter which 
determines the ratio between the coefficient of IZosiQ"^) and TZdQ^) , using the sum rule 
obtained by taking the derivative with respect to the Borel mass M^. We also eliminate the 
unknown coupling constant fy by taking the ratio between the obtained sum rule and the 
QSR expression for the vector current correlators in the vacuum. We thus get an expression 
for the spin-averaged scattering length ay as 

_ ttM^ r(3/{aM^) +t-f/{aM^) 

""^'^ ~ SmlJmN + mp,^) (1 + ^)(1 - e'^o/M^) + ^/jvf^ - c/M^' ^ ' 



with 



^ V " m'^ -L- ^ A^+'^ 



2 ^ m/, 

— rrin H 

27 ° 2 

t = ^{{uu){uu)N+{dd){dd)^)--m%Al-^^ + ^[B^^ + -B, + -B, 

— ix'^ ex 

b = 4n^mg{uu + dd) ^ (— C^), 

3 7r 



{uu) , 
o± 

where (■) denotes the vacuum condensate and 5*0 is the continuum threshold in the vacuum 
sum rule. The factors a, (3 and 7 appeared through the process of eliminating the parameter 
which determines the ratio between the residues of TZosiQ^) and TZdQ'^), and they are defined 
as a = 1 - e-(^o-^)/M^(l + /? = g + ^e'^o/M^ - Ae-(^o-™^)/M^, ^ = 1 + 

^ — (1 + -^)e~('^o-™^)A^^ with m = mp^uj- If we ignore TZdQ^) from the beginning, the 
corresponding formula is obtained by the replacement; a — > 1, /5 ^ 1 — e"™^/^^, 7 — > 1. 
In eq. (|^), we have used the following relations for the matrix elements as has been used in 
the study of QCD sum rules in the nuclear matter @, |10|: (i) ti'N cr-term Sj^jv is introduced 
through the relation mq{uu + dd) ^ = 2m7vE^jv (h) The nucleon mass in the chiral limit, mo, 
is introduced in favor of {^G'^)^ through the QCD trace anomaly: {^G'^)jq = — (16/9)mQ. 
(iii) Factorization is assumed for the vacuum four-quark condensates, (Qs") and (Q"*"), as is 
usually adopted in the literature 0. (iv) Factorization is also employed to estimate the 
nucleon matrix elements of the scalar four-quark operators {Q^)n and {Q'^)n after making 
the Fierz transform [^, i.e. {{qTXq)'^)^ ((^Q')^)Af — {qq) (qq) n ■ 

By repeating the same steps as above for J^, one gets the spin-averaged (f) — N scattering 
length as 

^ vrAf^ r,/3/(aM^)+t,7/(aM^) 

"'^ 3ml{mN + m^){l + ^){l-e-So/M^)-Qml/M^ + bs/M^-Cs/M^' ^' 
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with 

2 

Ts = ms{ss)N - + m%A'2, 

2 

hs = 87rV(ss) + ^(^G'2), 

o vr 

4487r^as , , 2 

Cs = {ss) , 

81 ^ ^ ' 

where the strange twist-4 matrix elements [i = 1 — 3) are defined similarly to the case of 
the p and u mesons. For the vacuum condensates and the quark masses in eqs. (|^) and (H), 
we use the standard values at the renormalization scale yU = 1 GeV 0: = 0.36, nig = 7 
MeV, rus = 110 MeV, {uu) = {dd) = (-0.28 GeV)^, (ss) = 0.8{uu). With these vacuum 
condensates and the continuum threshold Sq = 1.75 GeV^ for p, uj and 5*0 = 2.0 GeV^ for 
(j), the experimental values for rup^i^^^ are well reproduced. We thus fixed Sq at these values 
and use mp^^^ = 770 MeV and m<^ = 1020 MeV in eqs. (^ and (|). As a measure of the 
strangeness content of the nucleon, we introduce the parameter y = 2{ss) jy / {{uu) n + {dd) jv) 
and write {ss)^ = ymjqY,T,j^ /niq. For the nucleon matrix elements we use S^jv = 45 ± 7 
MeV, y = 0.2 and mo = 830 MeV obtained by the chiral perturbation theory [jl^]. Since 



we ignored the twist-2 gluon operators in eqs. (0) and (|), we consistently use the leading 
order (LO) parton distribution functions of Gliick, Reya and Vogt [jl2| to determine A""*""^ 
and {i = 2,4). It gives A^+'^ = 0.90, Al+'^ = 0.12, = 0.05 and Al = 0.002 at p^ = 1 
GeV^. For the twist-4 matrix elements Bi and B'^, we use our recent result extracted 



from the newest DIS data at CERN and SLAG. It is based on the SU(2) fiavor symmetry 
(i.e. B^ = {i = 1 — 3)) and a mild assumption on the matrix elements invoked by the 
flavor structure of the twist-4 operators. Both for the proton and the neutron, it gives 
Bi^ + B2/A + 553/8 = -0.24 ± 0.15 (-0.41 ± 0.23) GeV^ for the p {u) meson at p^ = 5 
GeV^ [Q. [Note that our p^ — N {N can be either proton or neutron) scattering length 
correspond to the isospin-spin averaged one.] 

Using these numbers for the matrix elements, the Borel curves for the p, (-j, 0-nucleon 
scattering lengths ap^^^fj, (eqs. (|^ and (||)) are shown in Fig. 1. We determined the values of 
Sq in order to minimize the slope of the curves at 0.8 < < 1.3 GeV^. They are 3.32 
GeV^ for p, 3.29 GeV^ for u and 4.40 GeV^ for 0. With the above parameters, r in eq. (|^) 
reads r = 0.04 — 0.05 -|- 0.40 = 0.39 GeV^ from the first to the third terms. Thus r is totally 
dominated by the twist-2 nucleon matrix element A'^'^'^ and the cancelling contribution from 
the first and the second terms makes the ambiguity in E^tv and mg less important. The 
t-term in eq. (^ reads t = 0.42 - 0.08 - 0.11 ± 0.07 (-0.18 ± 0.10) = 0.23 ± 0.07 (0.16 ± 0.10) 
GeV^ for p {u), which shows the contribution from the twist-4 matrix elements is sizable. 
To get an insight on the sensitivity of the results to the variation of t, we also showed 
0'P,uj without the twist-4 matrix elements in t with Sq = 3.35 GeV^. One sees that the 
inclusion of B^ reduces the ap^^ by about 20 % (30 %) for p {u). With the uncertainty 
in Bi in mind, we assign error bars as = 0.14 ± 0.07 fm and = 0.11 ± 0.06 fm, 
taking the values for cip,(^,</. around = 1 GeV^. For the case of a^, the value of ms{ss)N 
governs the whole result because of large m^, i.e. = 0.13 — 0.05 -|- 0.04 = 0.12 GeV^ and 
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ts = 0.066 - 0.003 + (twist -4 = 0) = 0.063 GeV^ from the first to the third terms in and 
tg. Due to the uncertainty in ms{ss)i\i, we read from Fig. 10,^ = 0.035 ± 0.020 fm. Some 
phenomenological analyses on the nucleon form factor and the nuclear force [|1^ suggest 
quite a large OZI violating (j)NN coupling constant g,j,NN / QuNN ~ 0.4. Equation (|) supplies 
a neat expression for the (pN (pN interaction strength in terms of the strangeness content 
of the nucleon, showing the importance of ms{ss)N rather than {s'^^j^Dys)]^- 

If we calculate the scattering lengths without TZdQ^) in eq. (|^), we get even smaller 
numbers for the scattering lengths: ~ 0.1 fm, ~ 0.08 fm and ~ 0.01 fm around 
M = 1 GeV. This way, the actual numbers for ap,^;,^ depends on the assumption made in 
the spectral function, although their typical order of magnitude does not change. 

One might be surprised by the smallness of these scattering lengths compared with a 
typical hadronic size (~ 1 fm). From eqs. (|^ and (H), one sees ay ~ l/{mN + my), since 
r and are dominated by the third and the first terms, respectively. If one applies the 
present method to the axial vector correlator, one can easily get the pion-nucleon scattering 
length in the isospin symmetric channel as a-^N- oc mNT^T,^ / {flijrijq + m^)), which is the 
same result as that of the current algebra Jl^. [In the chiral limit, o^at = 0, since S^tv = 0.] 
Therefore it is interesting to observe that our method of deriving the vector meson-nucleon 
scattering length is a generalization of the current algebra technique for the pion-nucleon 
scattering length. For the vector meson case, the common factor l/(m7v + fny) makes ay 
small. We believe this smallness of the V — N scattering lengths somehow sketches the real 
situation, although the actual numbers for ay are not trustable because of the simplified 
form for the spectral function in our calculation as was noted before. A model calculation of 
the p — N scattering amplitude based on an effective hadronic lagrangian suggests a similar 
small number for ap [jl8| . 



Let us finally discuss the mass shift of the vector mesons in the nuclear medium using the 
result for the scattering lengths here. In the dilute nucleon gas approximation, the \^-current 
correlator in the nuclear medium can be written as 

K^''i^,q)=^J rf'^e^^-(T(j;(x)jJ(0)))+E T^^?^^^'^^^' 



(27r)32/' 

By ignoring the Fermi motion of the nucleon, q = 0) can be approximated near the 



pol. ' 

J N.M. 

L'Clllil lliULlUli Ui LliC liUUlCUli, 1 

pole position as 




(10) 

where Amy = 12TTavPN{f^N + my)/miq with the nucleon density p^v From this relation. 
Amy can be viewed as a shift of m^ in the nuclear medium |jl9|. Our values for ap^u),<t> suggest 
that the effective mass for the vector mesons increases by about 27 — 57 MeV for p, 20 — 48 
MeV for uo and 5 — 13 MeV for at the nuclear matter density = 0.17 fm~^ ^0[. [Note 
that the validity of the mass shift discussed here hinges on the assumption that the off-shell 
energy dependence and the momentum dependence of the V — N scattering amplitude is 
weak within the range of the nucleon's Fermi momentum.] 
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The authors of ref. [1^] apphed the QSR method to study mass shifts of the p, u and 



mesons in the nuclear medium. Ahhough their approximation in the OPE side of the 
correlation functions is essentially the same as ours, eq. (|]), they found a serious decrease 
of these vector meson masses in the nuclear matter. Here we clarify the origin of this 
discrepancy. In the recent literature of the QSR method in the nuclear medium for baryons 



and mesons P, 10 1, the common starting point is that the density dependence of correlation 



functions is ascribed to the density dependent condensates: 

Il''-''-{q,PN) =Y.C^{q,p)m^i))p,, (11) 

i 

where Ci{q,p) and are the Wilson coefficient and a local operator, respectively, and 

we suppressed all the spinor and Lorentz indices. In the dilute nuclear medium, (Cj(/i))pj^ 
has been approximated as 

= (0.(/i)) + /^(O.(/x))^ + o(p^). (12) 



As is easily seen by inserting eq. (0) into eq. (|TT]), the approximation for the condensate 
eq. (0) is equivalent to the approximation eq. (P) for the correlator. Therefore the density 
dependent part of the correlator has to be analyzed from the point of view of the forward 
current-nucleon scattering amplitude as was done in this paper. In this approximation, what 
is relevant for the mass shift is the double pole structure at the pole position appearing in 
the forward amplitude. 



To understand the difference between our result and the one in it is convenient 
to recall the QSR for the vector meson in the vacuum. In the vacuum, the vector current 
correlator can be written as 11^!/ (g) = {qfj.qu — gfiuq'^)^{,q^)- As long as one uses a spectral 
function with a single narrow resonance and the continuum, both n(g^) and g^n(g^) can be 
used for the QSR analysis. The formula for the vector meson mass in the Borel sum rule 
obtained by using these two sum rules are, respectively 



rr^ _ (1 + ^)[1 - (1 + So/M^)e-^^/^'"] - DJM* + D^/M^ 

M2 ~ (1 + ^^)[1 -e-^o/^'] + -Z}6/2M6 

m|. _ 2(1 + f )[1 - (1 + Sq/M^ + 5'oV2M^)e-^°/^^'] - D^/M^ 

JP ~ (1 + ^^)[1 - (1 + S'o/M2)e-^oA^'] - DJM^ + D^/M^ 



(13) 



(14) 



where and Dq are the sums of the dim. =4 and dim. =6 condensates, respectively. Both 
sum rules give right amount for my. However, one should note that the sign of the power 
correction due to the condensates is opposite between the two sum rules. Thus the shift of 
the condensates due to the medium effects (the second term of eq. (^ or eq. ([T^)) is expected 
to cause opposite physical effects depending on which sum rule one uses. The authors of 
|T0| analyzed H^*^'(ci;, 0) = H^'^^''^(c(j, 0)/(— Scj^), with a simple pole ansatz for the vector 
meson (together with a scattering term) in the spectral function, which picks up the same 
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effect of the shift of the condensate as eq. ([T3|). On the other hand, if we recognize that the 
second term in eq. (|) is associated with the V — N forward amphtude through eq. (|]), we can 
easily see that it is uj'^Il^''^^'{u, 0) which has to be analyzed with a simple pole ansatz in the 
order 0{pn) as is shown in eq. ([10|) . In this case, the vector meson mass receives the effect 
of the change of the condensate as is expected from the formula eq. (|1^. In n^ *^'(c(j, 0) 



the density dependent part appears as a form of {pN/2mN)T{uj, 0)/a;^, which brings a factor 
-Amy{my/Q'^) (Q^ = -u^ > 0) instead of Amy in the first line of eq. (|10|). In this case, 
due to the additional factor l/Q"^, the double pole term can not be incorporated into the 
mass shift. Thus the use of eq. (|13D in the nuclear medium is simply wrong! Therefore an 



inadequate form of the spectral function in |]T0[ led to a fictitious "negative" mass shift. 



Since the second term in eq. has an unique relation with the — T-matrix around the 
pole position as is shown in eq. we believe the mass shift of the vector mesons in the 
nuclear medium in the context of the QSR should be understood as presented in this work. 

Finally, we wish to make a comment on the speculation on the in-medium behavior 
of the hadron masses existing in the literature. From the finite energy sum rule analysis, 
one gets the p meson mass as nip oc Hi^ip)]^^^ in the vacuum. Since K'^V')! decreases in 



the nuclear medium according to the formula eq. (|T2D, one might naively expect mp would 



also decrease [0]. We have illustrated, however, that a consistent organization of the QCD 
sum rule does not predict such a behavior. It would rather support (within our crude 
approximation) another naive expectation that a tendency of p — Ai degeneracy might occur 
in the nuclear medium, since the application of our method to the Ai meson gives decreasing 
rriAi- We also remind the readers that (1) |('^V)t| decreases as the temperature (T) goes up. 



while all hadron masses stay constant in the order OiT ) (2) a consistent organization 



of the sum rule at finite temperature certainly gives the same behavior |19] unlike the above 
naive expectation, and (3) this is because the pion-hadron scattering length is zero in the 
chiral limit. 

In conclusion, we have derived the p, u, (p -nucleon spin-isospin averaged scattering 
lengths ap,aj,(/) from QCD sum rules. We obtained very small positive numbers (corresponding 
to the repulsive interaction) for Op^^^^, ctp,,^,,/, ~ {rn^ + rrip^^^^), although the actual numbers 
depend on the factorization assumption for the nucleon matrix element of the scalar four- 
quark operator as well as the simplified form for the spectral function. This result suggests 
slight increase of these vector mesons in the nuclear medium, which is contradictory to the 
previous result by Hatsuda and Lee [|1^] . We have clarified the origin of this discrepancy 



and pointed out the problem of the analysis in ref . |I0[ . The detail of the calculation will be 
published elsewhere. 
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Figure Captions 



1 The Borel curves for the p, ui, 0-nucleon scattering lengths. The dashed hne denotes 
the one for p and u without the twist-4 matrix elements in eq. (^. 
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